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ABSTRACT 

We present kinematic measurements of thin and thick disk components in a sample of nine edge-on 
galaxies. We extract stellar and ionized gas rotation curves at and above the galaxies' midplanes using 
the Ca II triplet absorption features and Ha emission lines measured with the CMOS spectrographs 
on Gemini North and South. For the higher mass galaxies in the sample, we fail to detect differences 
between the thin and thick disk kinematics. In the lower mass galaxies, there is a wide range of 
thick disk behavior including thick disks with substantial lag and one counter-rotating thick disk. 
We compare our rotation curves with expectations from thick disk formation models and conclude 
that the wide variety of thick disk kinematics favors a formation scenario where thick disk stars are 
accreted or formed during merger events as opposed to models that form thick disks through gradual 
thin disk heating. 

Subject headings: galaxies: kinematics and dynamics — galaxies: formation — galaxies: structure 



1. INTRODUCTION 

The detailed distribution of stars in galaxies gives vital 
information regarding their formation and subsequent 
evolution. Of particular interest are the oldest stellar 
populations, which in the Milky Way are the thick disk 
and halo. These old components provide the best record 
of early gal axy assembly. Originally detect ed in edge-on 
SO galaxies (|Bursteinlll9"7l lTsikoudilll979( l. thick stellar 
disks have now been found in a wide variety of galaxies- 
SO's (Ide Griis &: van der K7^[1996I: Ide Griis fc Peletie3 
1997"; 'Pohlcn ct al." '2004), Sb's ("van der Kruit 1984; 

Shaw fc Gilmor c 198 ^; Ivan D o kkum et al. 199J; 

Morrison et"all 119971: IWu et all I2002D. and later 
type galaxies (|Dalc anton fc Bernstein! 120021: lAbe et al.l 



119991: iNeeser et al.l 2002: Yoachim fc Dalcantod l2006f l. 
Observations with HST have allowed thick disks 
in ot her galaxies to b e stud ie d as resolved p opula- 
tion (iSetli et all l2005l 120 07: Ti khonoy et"all 120051 : 
iTikhonov fc GalazutdinOTari2005; .Mould \2mm . while 
observations at high redshift s how potential thick disks 
in th e process of forming (|Elmegreen fc Elmegreenl 
l200^ . 

The most detailed studies of thick disks come from 
observations within the Milky Way. Since its discovery 
(|Gilmore fc Reidlll983D . the MW thick disk has been 
found to be structurally, chemically, and kinematically 
distinct from the thin disk. Structurally, star counts 
with large surveys such as SDSS and 2MASS reve al the 
galaxy is best fit with two disk components (e.g.. lOjhal 
120011: iJuric et al.l l2008f l. Chemically, thick disk stars 
are more met al-poor an d older tha n stars in the thi n 
disk (e.g., Rei d fc Maiew ski 1993; C hiba fc Beerg|l2000D . 
They are also significantly enhanced in a-elements, 
compar ed to thin disk stars of compar able iron abun- 
dance (iProchaska et al.l 120001: iT autvaisi ene et all 120011: 
iBensbv et al.ll2003l : iFeltzing et al.ii2003t iMishenina et all 
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200l iBrewer fc CarnevI [200l IBensbv et all [20051 : 
Brewer fc CarnevI l2006t iRami'rez et al.l l2007t) . Kine- 



matically, thick disk stars have both a larger velocity 

dispersion and slower net r ot ation than st ars in 

the thin disk (iNisse^ Il995l: IChiba fc Beerj l2000l: 
iGilmore et"all 120021: ISoubiran et all 120031 : iParker et all 
_2004t; iGirard et al"]l2006[ ). All of these facts lead to the 
conclusion that the thick disk is a relic of the young 
Galaxy. As such, it provides an excellent p robe of mod- 
els of d isk galaxy formation (see revi e ws by iNissen et al.l 
(|2003f l: lEYeeman fc Bland-Hawthorn! (|2002f )). 

Given these systematic differences between their prop- 
erties, thick and thin disks are likely to have dis- 
tinct formation mechanisms. The structure, dynam- 
ics, and chemical abundance of the thin disk strongly 
suggest that the majority of its stars formed gradu- 
ally from a thin rotating disk of high angular momen- 
tum gas (iFall fc Efetathioul 119801: [Chiappini et~al]ll997l: 
ICescutti et al.![2007!) . In contrast, the formation of the 
thick disk is still poorly constrained and is likely to be 
more complex. 

Thick disk formation models can be grouped into 
three broad categories. In the first, a previously 
thin disk is kinematic heated. In this scenario, stars 
form in a thin disk and increase their velocity dis- 
persion with time. This vertical heat ing can be 
rapid , due to inter a ctions and mergers (jQuinn et aLJ 
1993): 'Walker et al .l [199^ 'Velazq uez fc White! 119991: 



Chen e t al. 2001; Robin et al. 1996,) or gradual, due 
to scattering off giant molecular clo uds, spir al arms, 
and/or dark matter substructure (IVillumsen .19851: 
Carlberd [l987l: iHanninen fc Flvnnl 120021: iBensqn et al.! 



2004 iHavashi fc Chibal 120061 : iKazantzidis et al.! I2'007ir 
In the second formation scenario, stars "form thick" 
with star f ormation occurring above the midplane of 
the galaxy ([Brook et al.! |2004 ) or form with l arge ini- 
tial v elocity dispersions in large stellar clusters (jKroupal 
12002! ) ■ In the final class of models, thick disk stars 
are directly accreted from satellite galaxies. Nu- 
merical simulations have shown that stars in dis- 
rupted satellite galaxies can be deposited onto thick 
disk like orbits ([Abadi et al.! 12003! : [Martin et all [20041 : 
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Bekki fc Ch iball2001l: iGilmore et cd1[200l : iNavarro et all 
2004t IStat lcr 1988), pro ducing extended stellar debris 



such as seen around M31 (llbata et al.|[2005l : iKalirai et al.1 
l2006t iFerguson et all l2002f) . While these models were 
originally developed to explain the origin of the MW 
thick disk, they should work equally well for thick disks 
in other galaxies. 

Measuring the kinematics of thick disk stars is one of 
the best discriminators between the formation models. 
If the thick disk forms from a heated thin disk, we ex- 
pect the kinematics of the two components to be closely 
related. On the other hand, if the thick disk stars form 
outside the galaxy and are later accreted, we could find 
systems where the thick disk kinematics are completely 
decoupled from the thin disk. 

In this paper, we present observations of stellar and 
gas kinematics in nine edge-on systems as part of our 
continuing analysis of thick disks in a large sample of 
edge-on galaxies JD alcanton & Bernstein 2003). Com- 
pared to lYoachim fc Dalca ntonI (|2005l) . which presented 
the first two galaxies in this study, we have improved 
the analysis techniques and significantly expanded our 
sample size. 

2. OBSERVATIONS 

2.1. Target Selection 

We have carried out long-slit spectroscopic observa- 
tions using the Gemini North and South telescopes of 
nine galaxies drawn from the iDalcanton fc BernsteinI 
()200C0 sample of edge-on late-type galaxies. The original 
sample of 49 galaxi es was selected f rom the Flat Galaxy 
Catalog (Karachcn tsev et al.l Il993f) and imaged in B, 
R, and K^ (|Dalcanton fc BernsteinI [2000f) . This sample 
was selected to contain undistu rbed pure disk systems 
spanning a large range of mass. IDalcanton fc BernsteinI 
|2002!) used this imaging to demonstrate the ubiq- 
uity of thick disks around late- type galaxies, while 
lYoachim fc DalcantonI (j2006t ) used two-dimensional pho- 
tometric decompositions to measure the structural pa- 
rameters for the thick and thin disks. All the galaxies 
in the sample presented here have prominent thin star 
forming disks. 

Our spectroscopic program targeted galaxies spanning 
a wide range of masses (50 < Vc < 150 km s~^). The 
sample targets were limited to those that had thick disks 
that we believed we could isolate adequately~i.e., those 
that had significantly larger scale heights from the thin 
disk and that were bright enough that we could acquire 
spectra in reasonable observing times. This constraint 
caused several of the higher mass galaxies to be rejected 
from the kinematic sample, as the regions where the 
thick disk could be expected to dominate were simply 
too faint. This bias is consis tent with the conclusion of 
lYoachim fc DalcantonI (|2006f ) that the thick disk is more 
prevalent in lower mass galaxies. Our selection criterion 
limited the sample to ~20 galaxies of the original 49. 
We also selected galaxies to be at rcdshifts such that the 
Ca features did not land on night sky emission lines. In 
our initial observations, we submitted more galaxies than 
we could observe and let the Gemini observing special- 
ists select which galaxies would best fit with the queue 
scheduling. For the final observing runs we explicitly se- 
lected galaxies to ensure that a reasonable mass range 
was observed in the final sample. The properties of the 



final sample are listed in Table [T] 

2.2. Observing Strategy 

Based on the thin and th ick disk decompositions in 
'Yoac him fc DalcantonI ([2006") , we targeted regions of the 
galaxies where the flux is dominated by either the thin 
or thick disk stars. The two high est mass galaxies in ou r 
sample have notable dustlanes (JDalcanton et al.ll2004( l. 
and for these we offset the spectra slit to observe re- 
gions of the galaxy which should be optically transpar- 
ent. We discuss possible residual dust effects in detail 
in § |7| When selecting slit placement for the offplane, 
the direction of offset was based primarily on avoiding 
foreground objects and the ability to use a single guide 
star for all dither positions. 

For our instrumental setup, we used GMOS on Gem- 
ini North in longslit mode with a 0.5" slit and the 
R400_G5305 grating set to a central wavelength of '^ 
8440 A along with the OG515_G0306 filter. Similarly for 
observations from Gemini South, we used a 0.5" slit the 
R400+_G5325 grating and OG515_G0330 filter. For both 
GMOS setups, we binned the CCDs by 2 in the spatial di- 
rection during readout giving a pixel scale of 0.145"/pix 
in the spatial direction and 0.69 A/pixel in the spectral 
direction. The resulting spectra cover the wavelength 
range of ~ 6330 — 10570 A, although there is heavy resid- 
ual fringing redward of 9300 A. Exposure times for in- 
dividual frames were 900, 1200, or 1800 seconds. The 
midplanes were observed 3-5 times while offplane posi- 
tions were observed 18-51 times depending on the galaxy. 
Exposures were spatially dithered ^--^30" along the slit. 
These configurations allow us to simultaneously observe 
the Ha emission and Ca ii triplet absorption features out 
to large radii. 

All of the observations were executed in queue mode 
over five semesters. The observation details for each 
galaxy are listed in Table JH with details of the slit posi- 
tions listed in Table |3l 

2.3. Data Reduction 

A combination of Gemini IRAF packages, standard 
IRAF packages, and custom IDL code were used to re- 
duce our data. These procedures have been improved 
since initial results for EGG 227 and F GC 1415 were 
published in lYoachim fc DalcantonI (|2005( l and have been 
apphed to the entire data set. We bias corrected the im- 
ages using a fit determined from the overscan region fol- 
lowed by subtracting residual structure measured from 
a bias frame. Because both GMOS North and South 
are extremely stable, we were able to create average 
bias images by combining ~60 bias frames per observ- 
ing semester. We interpolated the three GMOS chips 
into a single image using the Gemini IRAF tasks, after 
which the standard IRAF reduction tools were used. For 
Gemini-South observations, we also needed to subtract 
a dark current correction of ~6-12 counts from the sci- 
ence frames. Gemini-North images showed no detectable 
dark current. Images were flat-fielded using GCAL lamp 
fiats that were taken every hour interspersed with the 
science observations, minimizing the amount of fringing 
present in the final frames. We applied a slit illumination 
correction using twilight sky observations. 

For w avelength calibrat i on, we used the night-sky at- 
lases of lOsterbrock et al.l ()1996D and lOsterbrock et al.l 
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TABLE 1 
Properties of Targeted Galaxies 





Galaxy Dist^ Vc 


It-R ^O.thin 


^0, thick ^thick/Lthin 




] 


FGC Mpc km s-i 


II 


II 


It 






227 89.4 106.0 


10.2 


1.8 


3.9 0.47 






780 34.4 75.0 


15.1 


3.1 


8.4 0.93 






1415 38.3 86.5 


18.3 


2.8 


6.6 0.95 






1440 70.9 150.5 


15.9 


2.3 


5.0 0.38 






1642 36.6 55.0 


12.5 


3.1 


10.0 0.19 






1948 36.9 54.5 


12.3 


1.6 


3.6 3.56 






2558 73.8 89.0 


9.2 


2.6 


3.6 0.47 






E1371 82.6 131.0 


7.7 


1.6 


3.4 0.27 






E1498 135.5 133.0 


7.6 


1.2 


3.8 0.19 




iRarachentsev et al. f2000) 












TABLE 2 
Observing Details 








Galaxy 


Gemini ID 


Observation Dates Midplane Exposure 


Offplane Exposure 


FGC 




# 


X time (s) 


# X time (s) 


1415 


GN-2003A-Q-6 


3-28-2003 to 06-06-2003 


3x900 




41x1200 


227 


GN-2003B-Q-51 


9-21-2003 to 11-22-2003 


3x1200 




27x1200 


1642 


GN-2004A-Q-54 


02-16-2004 to 06-24-2004 


3x1200 




51x1200 


780 


GN-2004A-Q-54 


02-20-2004 to 04-27-2004 


5x1200 




31x1200 


2558 


GN-2004B-Q-29 


07-15-2004 to 11-20-2004 


3x1200 




36x1200 


E1498 


GS-2004B-Q-44 


03-11-2005 to 06-10-2005 


3x1200 




50x1200 


1948 


GN-2005A-Q-21 


08-12-2004 to 08-24-2004 


5x1800 




18x1800 


E1371 


GS-2005A-Q-17 


04-05-2005 to 04-14-2005 


3x1200 




21x1800 


1440 


GS-2005A-Q-17 


02-11-2005 to 04-05-2005 


3x1200 




30x1800 






TABLE 3 












Slit Placement 










Galaxy 


Midplane Offset^ 


Offplane 


Offset 






FGC 


arcscc kpc 


arcsec 


kpc 


2/^0, thin '^ho,thick 


227 




0.0 0.0 


3.0 


1.3 


1.7 0.8 


780 




0.0 0.0 


6.5 


1.1 


2.1 0.8 


1415 




0.0 0.0 


5.4 


1.0 


1.9 0.8 


1440 




0.5 0.2 


4.5 


1.5 


2.0 0.9 


1642 




0.0 0.0 


4.2 


0.7 


1.4 0.4 


1948 




0.0 0.0 


3.1 


0.6 


1.9 0.9 


2558 




0.0 0.0 


3.9 


1.4 


1.5 1.1 


E1371 




0.5 0.2 


2.8 


1.1 


1.8 0.8 


E1498 




0.0 0.0 


2.0 


1.3 


1.7 0.5 


^Midplane offset to avoid obvious dust lanes. 









()1997f ) to create a sky hne hst containing only hnes (or 
stable unresolved doublets) that could be centroided with 
our instrumental set-up. For each science exposure, we 
identified 100-flO sky lines to use for rectification. We 
then used these lines for a 5th order Legendre polynomial 
fit for wavelength calibration, and rebinned our spectra 
to a common dispersion. Typical dispersions were 0.69 
A pixel" ^ with calibration arc lamps showing a FWHM 
of 3.8 A. The wavelength solutions were stable over each 
observing night. 

Sky subtraction proved difficult because of the large 
number of strong sky emission lines. If we use stan- 
dard sky subtraction techniques, we find that there 
are large systematic residuals left on our frames due 
to variation in the width of the slit along its length. 
The RMS deviation in the centroid position of a sin- 



gle sky hne is --0.07 A while the RMS of its Gaus- 
sian FWHM is 0.11 A. This is a surprisingly high vari- 
ation for the width of the slit. We have tried the sky - 
subtraction techniques described in lKelson et al.l (|200Gl ) 
and find that the system atic residuals remain, although 
the iKelson et al.l (|200GD sky-subtraction technique does 
eliminate problems associated with wavelength rectifi- 
cation and interpolation. Having eliminated our data- 
reduction procedure as the cause, we conclude the high 
dispersion in sky line FWHM is indicative of a system- 
atically varying slit width. In many cases, such residu- 
als can be removed using the nod-and -shuffle technique 
(jGlazebro ok fc Bland-HawthornI [20011 ). Unfortunately, 
our galaxies are too large (~1 arc minute, or 1/3 of the 
total slit width) to make effective use of traditional nod- 
and-shuffic. 



Yoachim & Dalcanton 



To remove the systematic residuals present in the 
bright sky Unes, we employ a nod-and-shuffle like tem- 
plate subtraction. Because we placed different galaxies 
on different spatial sections of the chips, all of the slit was 
illuminated by sky for at least some observations. We 
therefore could construct high S/N sky frames by mask- 
ing objects in our 2-d spectra and combining the wave- 
length rectified frames. By doing this, we create a deep 
sky frame for each observing quarter. We then remove 
the sky background by selecting a sky-dominated region 
in a science frame and scaling the sky image column- 
by-column to match the science frame sky region, then 
subtract the rescaled sky frame from the science im- 
age. In most cases, we were forced to apply sky frames 
generated from different observing semesters to the sci- 
ence frames. Luckily, our instrument setup quarter-to- 
quarter was identical, and the GMOS instruments are 
stable enough that this technique works well at remov- 
ing systematics caused by the variable slit width. This 
sky subtraction technique appears to give results compa- 
rable to nod-and-shuffle technique for individual frames. 
Our sky subtraction procedure incurs a small signal-to- 
noise penalty, but is effective at removing the systematic 
residuals from moderate sky lines (Figure [T]) . 

This excessive agonizing over sky subtraction is de- 
manded by the very low surface brightness levels of our 
targets. For an individual midplane image, the bright- 
est part of the galaxy is ~ 20% brighter than the sky 
level, and for individual offplane images the signal is only 
~ 11% the sky background. Examples of the spectra ex- 
tracted over the central 14" spatial extent of the galaxy 
before and after sky subtraction are shown in Figure [2l 

When we combine several hours of observations we 
are more sensitive to low surface brightness features, 
and find some wavelengths are still dominated by sys- 
tematic noise. Even with our sky template correction, 
some sky lines are so bright that some systematic resid- 
uals remain. When we use conventional sky-subtraction 
techniques, residual errors have maximum deviations of 
±55% while the template subtraction gives deviations of 
±38%. While deviations of 38% swamp out the signal 
from any stellar absorption lines near bright sky lines, 
the residual deviations for smaller sky lines are decreased 
to a level where the stellar absorption lines can be accu- 
rately measured. In Figure (TJ we compare the two sky 
subtraction routines. The extracted spectra look simi- 
lar, with both being dominated by the sky line residuals 
redward of 8750 A. The template subtraction is able to 
eliminate the residuals left from the sky line at 8555 A, 
just to the right of the weakest Ca ii triplet line, and 
reduces the large residuals at the reddest wavelengths 
plotted. 

After the sky had been removed, the images were 
Doppler-corrected for motion relative to the Local Stan- 
dard of Rest and combined. Before cross-correlation was 
performed, the spectra were rebinned into logarithmic 
wavelength bins. 

3. ROTATION CURVES 

3.1. Ha Rotation Curves 

Both our midplane and offplane observations show 
strong Ha emission. For each galaxy, we extracted a 
series of 1-D spectra by summing 28 pixels (~ 4") along 
the spatial dimension. The ionized gas rotation curve 



was fit with a Gaussian peak to the Ha line. In princi- 
ple, an envelope-tracing method would produce a more 
robust measure of the rotation curve. However, we find 
that the width of the Ha hues (FWHM~3.8 A) are iden- 
tical to the instrumental dispersion as measured from 
the arc lamps (FWHM'^3.8 A), and we would thus not 
gain much accuracy from a more detailed rotation curve 
extraction. 

The [Nil] and [Sii] lines are present as well, but the Ha 
line is so strong that we found no additional advantage 
in fitting all the emission lines simultaneously. We find 
typical uncertainties in the central wavelength of the Ha 
Gaussian peak of 1-2 km s~^ for midplane observations 
and 4-7 km s^^ for offplane observations. 

To double check the accuracy of our extracted rotation 
curve, we fit rotation curves to night sky lines before the 
background is subtracted off. Perfect calibration would 
result in sky line rotation curves with zero rotation. The 
central wavelengths of the sky lines vary with an RMS 
error of 2.4-3.5 km s~^, with the higher value resulting 
from larger spatial extraction windows. Most of this scat- 
ter can be attributed to uncertainties in the wavelength 
rectification solution. With fewer sky lines around Ha 
compared to the redder regions of our spectra, the rec- 
tification is not as well constrained. Overall, these tests 
suggest that we are able to extract the ionized gas rota- 
tion curve with an error of a few km s~^. 

The resulting Ha rotation curves are plotted as solid 
lines in Figure [H Our data show a tight agreement be- 
tween the midplane and offplane Ha curves, which is a 
good sign that dust is not obscuring the midplane rota- 
tion curves. If we were observing along major dustlanes, 
we could expect to see the offplane observations rotat- 
ing faster than the midplane, especially at small galactic 
radii (see §[7]). 

We leave a detailed analysis of the gas kinematics for 
a later paper. At this time, we simply note that the 
midplane and offplane Ha rotation curves are surpris- 
ingly well matched. This is slightly unexpected, as sev- 
eral recent studies have found extended gaseous halos 
of edge-on galaxies to be lagging i n rotational speed - 
when compared to the midplane gas (iHeald et al.|[2006al : 
iHeald et al.ll2007t iFraternah fc Binnevll2006[ ). These off- 
plane lags have been detected in both the diffuse ion- 
ized gas (DIG) and HI. There is some difficulty in com- 
paring our measurements of longslit rotation curves to 
other detailed measurements of offplane gas which typi- 
cally utilize 2-d information from radio (iBarbieri et al.l 
I2005t IFraterna h & Binncv 2006), Integral Field Units 
(jHeald et a i"2006a: Hcald et al.ii2007i) . and Fabry-Perot 
spectra (JHeald et al. 2006b) all of which detect gas at 
larger scale-heights than those probed with our offplane 
measurements. The other major difference between these 
previous studies and our offplane rotation curves is that 
we have targeted lower mass galaxies. The studies cited 
above target galaxies with 220 > Vmax > HO km s~^ 
while the sample studied here extends to galaxies with 
rotation speeds of less than 60 km s~^. 

The gaseous lags observed in other systems are usually 
modeled with either a galactic fountain that ejects gas 
to large scale-heights or with a gas infall model where 
galaxies slowly accrete rotating gas. The lack of signifi- 
cant lags in our Ha rotation curves could simply be a sign 
that these galaxies are not as active in forming galactic 
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Fig. 1. — Results from different sky subtraction techniques for the midplane of FGC 1415. On the left, we show the results from standard 
sky subtraction techniques and the right panels show our sky template subtraction. Top panels show the raw galaxy spectrum before 
the sky has been subtracted. Middle panels show a single subtracted frame and the final combined image. The bottom panel shows the 
combined spectrum summed along the spatial dimension. An arrow points out a sky line residual present in the standard subtraction that 
is eliminated in template subtraction. The brightest sky lines leave large residuals in both cases, but the magnitude of residuals is decreased 
significantly with the nod-and-shuffle-like technique (see the lines near ^ 8770 A, for example). 



fountains or accreting gas as the more massive galaxies. 
3.2. Ca II Rotation Curves 

To derive absorption line rotation curves, we require 
higher signal-to-noise than for the Ha rotation curve. 
We therefore sum the 2D spectra in the spatial direction 
until the ID spectra reaches an adequate S/N (~ 15 per 
spectral pixel). The resulting bins have variable widths 
across the face of the galaxy, but roughly comparable 
S/N per bin. For the central regions of the galaxies the 
bin size is around 10" while the outer regions and offplane 
components have bin sizes '^ 20". These bins correspond 
to ~3-6 kpc at the typical distances of the galaxies. For 
reference, the typical exponential disk radial scale lengths 
are h^ ~ 12". 

Extracting kinematic information from this 
data required develo p ing a new procedure. In 

lYoachim fc DalcantonI (|2005f ). we tried both direct 



X^-fitting of a template spectrum as well as cross- 
correlation of the galaxy with a stellar template to 
measure the stellar rotation and line-of-sight velocity 
dispersion (LOSVD). We have since concluded that 
these traditional methods are not optimal for our data. 
Direct fitting of a template star results in the template 
being over-broadened (i.e., the fitted LOSVD diverges 
to large values). This can be understood as the template 
star fitting the continuum region of the galaxy spectrum 
at the expense of a small portion of the absorption line. 
Because the normalized continuum is very low S/N, it 
is best fit by a straight line, which is equivalent to a 
stellar spectrum which has been smoothed b y a very 
broad filter. In lYoachim fc DalcantonI (|2005( ) we were 
forced to hold the velocity dispersion fixed during the x^ 
minimization to prevent this problem. Cross-correlation 
is also problematic, as the bright sky lines leave regions 
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Fig. 2. — Examples of spectra before and after sky subtraction. The top panel shows the results of a single midplane exposure before 
and after sky template subtraction (top and middle curves respectively). The middle panel shows a single ofEplane exposure before and 
after extraction. Dotted lines show the RMS noise level in the spectra. The bottom panel shows the final midplane and offplane spectra 
after all the frames have been averaged together. The largest systematic residuals from the sky lines have been masked. The three vertical 
marks show the location of the Ca triplet absorption lines. All of the spectra were extracted over the central 14" of the galaxy. 

trying different template stars and find our uncertainties 
are never dominated by template mismatch. 

Because we have modified the traditional cross- 
correlation technique, we have no formal means of cal- 
culating uncertainties in our fitted velocity and LOSVD. 
We therefore run a series of Monte Carlo realizations to 
quantify the errors in our fitting procedure. For each 
galaxy, we create 100 artificial 2D spectra. A template 
stellar spectrum is shifted to match a realistic rotation 
curve, and broadened to simulate both stellar velocity 
dispersion and instrumental resolution. We vary the de- 
tailed shape of the rotation curve and velocity dispersion 
for each realization by ^ 20%. The fake spectra have ra- 
dial exponential fiux profiles similar to the real galaxies. 
We add Poisson noise to the artificial spectra, as well as 
systematic residuals by adding regions of sky from our 
science frames that do not have any detectable objects. 
Thus, our artificial spectra have both the same Gaussian 
sky background and similar systematic residuals as the 
real data. 

Once the artificial spectra are made, we extract and an- 
alyze ID spectra identically to the real data (i.e., we use 
the same extraction windows and the cross-correlation 
with masking procedure). In many instances, we found 
that our measured LOSVD poorly matched the input. 
The loss of reasonable LOSVD measurements is domi- 
nated by how many of the Call lines are masked due to 
sky line contamination. We therefore clip points where 
the Monte Carlo error analysis suggests we cannot reli- 
ably recover the input parameters (i.e. the RMS error 



of very low S/N and systematic residuals caused by vari- 
ations in the slit- width (Figured]). Without a constant 
S/N throughout the spectra, the cross-correlation peak 
can become skewed by noisy regions. 

To extract both velocity and velocity dispersion in- 
formation from our spectra we developed a modified 
cross-correlation technique that allows regions of very 
low signal-to-noise to be masked. This modification pre- 
vents us from using the usual mathematical techniques 
involving Fourier transforms and instead utilizes a brute- 
force methodology. What it lacks in mathematical ele- 
gance, our procedure makes up for in functionality by be- 
ing the only procedure we know of that works on spectra 
that are both low S/N and contaminated with systematic 
residuals. We describe our modified cross-correlation in 
detail in Appendix [X] and compare its results to more 
traditional analysis methods in Figure [9l It may also 
be possible to use a penalized pixel-fitting technique to 
measure the kinematics from our spectra, but simulations 
show that the fitted parameters can become biased when 
the S/N is low (60), or the L OSVD is poorly sampled 
(jCappellari fc Emsellemll2004f ). 

For the stellar template, we used a Kill spectrum 
of star IID4388 downloaded from the Gemini archive 
along with accompanying calibration frames of program 
GN-2002B-Q-61. The stellar spectrum was reduced and 
extracted using the Gemini IRAF routines. Once ex- 
tracted, the ID stellar spectrum was broadened with a 
Gaussian kernel to match the instrumental resolution of 
our observations. We found no significant changes when 
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Fig. 3. — An example of our extracted galaxy spectra. The solid line shows the normalized galaxy spectrum. Red regions mark where 
the spectra was masked due to sky line contamination. The noise spectrum (multiplied by 5) is plotted as a dotted line. The blue dashed 
line shows the best fit shifted and broadened stellar spectrum. 

between input and output is > 50 km s^^ or the output 
has a systematic error of > 20 km s^^). These chpped 
regions typically correspond to regions of the rotation 
curve where the Ca triplet line passes through a large 
sky residual. 

Our final extracted rotation curves, LOSVDs, and 
Monte Carlo derived uncertainties are plotted in FigureH 
along with R-hand images of the galaxies showing the 
Gemini longslit placements. 



4. STELLAR KINEMATICS 

Although we attempted to place our slits in regions of 
the galaxies where the thin and thick disk light makes 
up the majority of the flux, it is nearly impossible to 
target regions where one stellar component completely 
dominates the flux. In the lower-mass galaxies, we found 
that the thick disk is a major stellar component and we 
should expect spectra taken along the midplane to in- 
clude a large amount of thick disk light. In the higher 
mass galaxies, the thin disk is the dominant component, 
and we are forced to observe off-plane regions that still 
contain a large fraction of thin disk light. Us ing the pho- 
tometric fits of lYoachim fc DalcantonI ()2006[ ). we can es- 
timate the fractional flux levels of the thin and thick disk 
at each slit position. Because each slit position should in- 
clude both thin and thick disk stars, we make an attempt 
to model the true underlying rotation curves for each 
population. For simplicity, we assume that the thin and 
thick disk stars are each rotating cylindrically and there- 
fore have the same rotation curve for both the on and 
off-plane observations. We discuss this choice in more 
detail in [J6] 

The details of the vertical profiles of the stellar disks 
(exponential vs sech^) can dramatically influence what 
fraction of the midplane light bel ongs to thin disk stars 
versus thin disk stars. As in lYoachim fc DalcantonI 
([2005'), we adopt a series of photometric decomposition 
models that should cover the full range of possible thin 
and thick disk fractions. At one extreme, we use a sim- 



ple model where we assume the midplane is composed of 
only thin disk light and the offplane observations purely 
thick disk stars. As a more accurate model, we use 
the thin/thick fraction s from the best fitting models of 
lYoachim fc DalcantonI ()2005l ) as well as models where 
we vary the parameters by their l-cr values to create a 
"bright-thick and faint-thin" model along with a "faint- 
thick and bright-thin" model. The differences between 
the thin and thick disk scale lengths are small enough 
that we do not expect much radial variation in the frac- 
tion of thin and thick disk l i ght. 

In lYoachim fc DalcantonI (J2005), we fit analytic func- 
tions to the stellar rotation curves to decompose the thin 
and thick disk components. This worked well for the ini- 
tial two galaxies we observed, but our expanded sample 
now includes galaxies with slightly irregular kinematics 
that are not well described by common parameterizations 
of rotation curves. Instead of using an analytic function, 
we use the midplane Ha rotation curve as a basis function 
for the overall shape of the rotation curve. Because we 
are most interested in finding the velocity of the thin and 
thick disk stars relative to each other, we compare them 
both to the well resolved and high signal-to-noise mid- 
plane Ha rotation curve. Using the Ha rotation curve 
reduces the number of parameters that need to be fit to 
characterize the stellar rotation curves. 

We model the stellar rotation curves as Vstars{R) = 
xVHa{R) + c. We constrain c to be in the range ±5 
(to account for any small error in wavelength calibration 
between frames or regions on the chip) and x is limited 
to —1 < a; < 1.4, allowing for stars to be rotating faster 
than the gas by up to 40% (a;=1.4), not rotating (x—O), 
or counter-rotating with the opposite velocity of the Ha 
(x=-l). 

The decomposed rotation curves are plotted in Fig- 
ure [5] The left hand panels show the best fit stellar 
rotation curve scaled from the Ha at each slit position. 
If there were no cross-contamination of thin and thick 
disk stars, then the offplane and midplane rotation curves 
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Fig. 4. — Rotation curve measurements for each galaxy. Top panels: /?-band images of each galaxy. The color scale goes from dark blue 
{flu = 21) to gieen(fiii = 23), to red/white {/ift = 25.5). Solid black lines have been drawn where the Gemini long-slit jaws were placed. 
Middle Panels: Rotation curves for midplane (blue) and offplane (red). Points with error bars are from Ca II measurements. Vertical error 
bars are uncertainties derived from Monte Carlo simulations, horizontal error bars show the spectral extraction regions. Small lines show 
velocities measured from the Ha emission lines. Bottom Panels: Stellar velocity dispersions measured from the Ca II feature. All error bars 
are from a Monte Carlo simulation. Points with overwhelmingly large error-bars or large systematic uncertainties have been omitted. 



would show the true thick and thin disk kinematics. The 
right hand panels show the more realistic case where we 
have adopted likely amounts of thin and thick disk con- 
tamination at each slit position before inferring the un- 
derlying kinematics of each population. 

For the higher mass galaxies, we find no substantial dif- 
ference between the thin and thick disk rotation curves, 
even when we correct for the expected cross contamina- 
tion. There is a slight tendency for the thick component 
to be lagging, but never by more than 5 km s~^. In the 
higher mass galaxies, we have therefore either failed to 
observe an offplane region with a high enough thick disk 
flux fraction, or the thick disks are not lagging signifi- 
cantly compared to the thin disk in these systems. 



For the low-mass galaxies, we find a wide range of be- 
havior. The fits for FGC 1948 diverge, as the stellar 
rotation curves do not show coherent rotation at either 
slit position. For the rest of the galaxies, the best fits 
find thick disks that are slightly lagging compared to the 
thin (FGC 2558, FGC 1415), that are lagging to the ex- 
tent of near non-rotation (FGC 1642, FGC 780), and 
that are fully counter-rotating (FGC 227). We note that 
th ere is strong qualitativ e agree ment with initial results 
in lYoachim fc DalcantonI ([2005h for FGC 1415 and FGC 
227. 

4.1. Velocity Dispersions 

The low signal-to-noise of our spectra prevents us from 
reliably measuring velocity dispersions for many of our 



The Kinematics of Thick Disks 



FGCE 1498 133 km s' 



FGCE 1371 131 km s' 



1 » 


~ -5 


9750 


-r 9700 


;: 9650 


J 9600 


■- 9550 


> 9500 


^fSlJ 


^ 80 


a 60 

A 40 


'° 20 










y 

s I? 
~ -5 =■. 







^ 6000 

I 5900 

— 5800 
> 

5700 
^ 150 
g 100 

e 50 





i' } 



-10 10 

R (arcsec) 



-10 10 

R (arcsec) 







FGC1440150 5kms 


1 




10 


V ' 


^:^,^j;-^v^vr^j. • ' 


-^ 


" 




—^^"""^"^^^^^■^ 


: 


^ -5 
-10 




'^^^ ^^^' 


] 


4400 
- 4350 
'" 4300 
E 4250 
w 4200 
> 4150 

4100 


t -^ 


^y^^^^ 


1 


100 

"«, 80 
^ 60 
« 40 
>= 20 



- 


H«H ' 


-! 



-60 -40 



-20 20 40 

R (arcsec) 



60 



Fig. 4. 



continued. 



galaxies. Most of the galaxies with high quality spec- 
tra have very low velocity dispersions, as we would ex- 
pect from systems predominantly supported by rotation. 
Given that our instrumental resolution is 60 km s~^, we 
are unlikely to resolve the line widths in galaxies where 
<j/Vc < 0.6, for the Vc < 100 km s^^ galaxies that dom- 
inate our sample. The major exceptions are FGC 1948, 
which has an irregular rotation curve, and FGC 227, 
which has a counter-rotating thick disk. 

FGC 1948 has surprisingly large LOSVD, with many 
regions of the disk having cr > 100 km s~^. For com- 
parison, most of the other galaxies in our sample have 
LOSVDs across the disk of ^ 50 km s~^, essentially the 
same as our instrumental resolution at the Ca ii triplet. 
The stellar rotation curve for FGC 1948 also shows large 
deviation from the Ha RC, suggesting that the stars 
in this galaxy might not be fully rotationally supported 
and/or fully dynamically relaxed. 

FGC 227's LOSVD also deviates from the simple in- 
terpretation of a dynamically cold rotating disk. In the 
midplane observations, the central regions of FGC 227 
appear cold {a ~ 40 km s~^), but the outer disk reaches 
LOSVD values of 100-150 km s~^. This makes little sense 
for a galaxy with a well defined rotation curve as the in- 
trinsic stellar velocity dispersions should be decreasing 
with radius. In contrast, the LOSVD can be well ex- 
plained by a rotationally supported galaxy if there are 
two stellar populations moving in opposite rotational di- 
rections. As our rotation curve decomposition showed, 
FGC 227 is best fit by a model where the thick disk is 
co unter-rotating relative t o the thin disk. As we showed 
in lYoachim &: DalcantonI ()2005D , this would cause an in- 



crease in the observed velocity dispersion of order 50 
km s~^. Similar projection effects are found in elliptical 
galaxies with counter-rotating cores as they also show 
radially increasing LOSVDs (jGeha et al.ll2005D . 

5. HOW MUCH COUNTER ROTATING MATERIAL COULD 
THERE BE? 

Inspired by the best-fit rotation curve for FGC 227, 
we investigate the possibility that all thick disks con- 
tain some fraction of counter-rotating stars. Our data is 
able to place tight constraints on the amount of counter- 
rotating material since both the offplane rotation curves 
and the midplane LOSVD will be strongly affected by 
any counter-rotating stars. 

In Section 31 we imposed thin and thick disk flux frac- 
tions based on previous photometric decompositions. We 
now leave the flux fractions as free parameters and in- 
stead hold the rotation curve shapes fixed. We fit two 
simple models, each with two kinematically independent 
stellar components. In the first model, we assume there 
are two stellar components, one rotating identically as 
the gas and one with zero net rotation, as one might ex- 
pect for a stellar halo. The final observed rotation curve 
is a flux weighted average of these two curves and we fit 
for the best fitting flux ratio. We restricted the explored 
parameter space such that the rotation curves had to 
be some positive linear combination of the midplane Ha 
and a non-rotating or counter-rotating rotation curve. In 
the second model, we assume the second component is 
counter-rotating with a velocity one-half the magnitude 
of the Ha rotation curve. For both models, we calculate 
uncertainties from the covariance matrix and scale them 
upwards such that the reduced-x^ equals unity (i.e., we 
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Fig. 5. — The results of fitting various rotation curve models to our data. Tlie top left panels show fits of the simple model where 
the midplane and offplane observations are fit independently. Upper right panels show shaded regions show the range of fits derived from 
varying the fraction of thin and t hick disk light at each slit po sition. Solid lines show the fits for when we use the thin and thick disk 
fractions of the photometric fits in I Yoachim fc DalcantonI (120061 ). Lower panels show the observations as points and solid lines show the 
models from the above panels once they have been flux weighted and binned in the same manner as the observations. Throughout, red 
is used for thick disk/offplane and blue is used for thin disk/midplane. Each panel has a dashed line showing the W50/2 value from the 
literature. FGC 1440 is not shown because we failed to measure a stellar rotation curve in the offplane position. 

assume our model should be a good fit). We do not cal- 
culate uncertainties when the fit converges to a boundary 
condition. We also do not construct detailed models for 
cases like FGC 1415 where the stars could be better fit 
with a faster rotation curve than the gas; these galaxies 
naturally converge on the boundary condition of hav- 
ing no second component. It should be emphasized that 
these are simple toy models, and we have no direct ev- 
idence of counter rotating thick disk stars beyond the 
strange rotation curve of FGC 270. For example, if we 
observed a MW hke (Vc — 220 km s~^) galaxy that had 



a 10% (by flux) thick disk lagging at 40 km s~^, we would 
compute a maximum counter-rotating fraction of 1% and 
a non-rotating fraction of 2%, despite all the stars being 
co-rotators. 

The resulting fractions of non-rotating of counter ro- 
tating stars are plotted in Figure [6] and are listed in 
Table 2] The midplane stellar rotation curves are typ- 
ically consistent with the Ha rotation curve, with 6 of 
the 9 midplanes being best fit without a non-rotating or 
counter-rotating component. The remaining three galax- 
ies do have midplane rotation curves that are consistent 
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TABLE 4 
Non-Rotating and Counter Rotating Fractions 



FGC 


Non-Rotating Fraction 


Counter-Rotating Fraction 




Thin Disk 


Thick Disk 


Thin Disk 


Thick Disk 


227 


O.lit 0.1 


0.8± 0.1 


0.1± 0.0 


0.5± 0.1 


780 


0.0 


0.3± 0.2 


0.0 


0.2± 0.3 


1415 


0.0 


0.0 


0.0 


0.0 


1440 


0.0 




0.0 




1642 


0.0 


0.4± 0.1 


0.0 


0.3± 0.1 


1948 


0.3± 0.2 


0.1± 0.2 


0.2± 0.1 


O.Oit 0.5 


2558 


0.2± 0.0 


0.4± 0.1 


0.2± 0.0 


0.3± 0.1 


E1371 


0.0 


0.0 


0.0 


0.0 


E1498 


0.0 


0.0 


0.0 


0.0 



with the presence of an additional lagging component. 
FGC 1948 is low mass with a surprisingly large LOSVD. 
FGC 227 is the counter rotator with a LOSVD that dra- 
matically increases with radius. FGC 2558 is the only 
galaxy to show a large discrepancy between midplane 
and offplane Ha rotation curves, has a stellar lag that 
appears to be only on the receding side of the galaxy. 

The offplane spectra show larger evidence for non- or 
counter-rotating motion, with only 3 of the 9 galaxies 
requiring no slow rotating component. This effect can 
be seen in Figure [SI where all of the offplane spectra 
show a preference for equal or larger value of the counter- 
rotating fraction than seen in the midplane. 

6. EXPECTED LAGS 

Having found a wide range of thick disk behaviors, we 
now investigate the expected stellar lags we should see 
in our sample of thick disks using a dynamical model 
originally designed for the Milky Way. The large scale 
height of thick disk stars implies they have larger veloc- 
ity dispersions than thin disk stars. If the larger vertical 
velocity dispersion also reflects a larger radial velocity 
dispersion, then the larger random motions of thick disk 
stars should lead to their requiring less rotational sup- 
port. The thick disk stars should therefore lag in velocity 
compared to the kinematically colder thin disk stars and 
io nized gas. 

iGirard et al.l ()2006[ ) use the Jeans equation and a series 
of reasonable assumptions to model the expected thick 
disk lag in the MW as a function of height above the 
midplane. While this model was built to explain the 
observed lag of thick disk stars in the Milky Way, the 
formalism is easily generalizable to the galaxies in our 
sample. 

Using the Jean's equation, iGirard et al.l ()2006l ) find 
that the rotational velocity of a thick disk rotating in a 
Plummer dark matter potential with an embedded thin 
disk is given by: 



+0.5Afl- 

(i?2 +Z^+ a2)3/2 



Ta,bR (1) 



a2)3/2 



1-^ 
4 



^disk 



{R,z), 



where R, z, and are galactocentric cylindrical coor- 
dinates. The term vq is the average thick disk velocity 
in the direction of galactic rotation, a^ and cfq are the 
radial and tangential components of velocity dispersion 



for the thick disk stars, Vc is the local standard of rest 
velocity at the radius of interest, Vdisk is the portion of 
the thick disk rotational velocity due to the gravitational 
potential of the thin disk, hz^^-^^. is the exponential thick 
disk scale height and a is the halo core radius. The term 
Ta,b lets one approximate the thick disk as entirely self 
gravitating, or gravitationally dominated by the embed- 
ded thin disk. Because the thick disk mass is small com- 
pared to the total gas and thin disk mass in all of our 
galaxies, we choose to use T;, ~ 2/hji. The A term takes 
values of 1 or in order to include or exclude the velocity 
dispersion cross-term. 

We calculate dynamical models for three fiducial 
galaxy masses and three thick disk velocity disper- 
sions. We use realisti c galactic parameters taken from 
lYoachim fc DalcantonI ((2006) to generate /izthicfc ^-^id 
^disfc(-^'O)- ^o^ terms for which we do not have ex- 
plicit measurements, we use t he approximation a[k pc] « 
13(hrt^,Jkpc]/5)i-°5 given bv iDonato eTall (|2004f ). as- 
sume (Te ~ (JR, and set A = 0. We compute models for 
different values of an, as this is the dominant term in pro- 
ducing stellar lags. For simplicity, we assume the thick 
disk velocity dispersion does not vary with height above 
the midplane. This l ast approximation is not particu- 
larly valid given that 'Girar d et al.l (|2006f ) find that the 
velocity dispersion in the MW increases with a slope of 9 
km s^^kpc^^. However the difference between a variable 
and constant velocity dispersion will be most pronounced 
at large scale heights, beyond the range probed by our 
observations (z ~ 1.5 — 2 kpc). The resulting models are 
plotted in Figure [7] along with the lags we have measured 
in our galaxies. For reference, we also include a model 
using the same assumptions but with morphology and 
velocities similar to the Milky Way in Figure [T] 

For most of the galaxies where we measure a thick 
disk lag. Figure [7] shows the thick disk kinematics could 
be well explained by a population with radial velocity 
dispersion of between 15 and 30 km s~^and Vc/cr < 4. 
As before, the major exception is FGC 227. The stellar 
lag for FGC 227 is so severe that it would imply the 
thick disk is completely supported by random motions. 
However, we only detect flattened stellar populations in 
FGC 227, again consistent with our interpretation that 
the thick disk is counter-rotating in this system. 

To verify that our model galaxies are reasonable, we 
use an identical procedure to build a MW-like model. 
Our MW-like model is a fair fit to actual observations of 
the MW. The measured thick disk velocity dispersion in 
the solar neighborhood is 50 km s~^, for which our model 
correctly predicts the midplane thick disk lag of 30 km 
s^^. On the other hand, the increase of the thick disk lag 
with scale height is poorly fit by our model; the observed 
lag increases with a slope of 30 km s^^kpc^^, and our 
model has a slope around half that. This is purely due 
to our choice to hold the velocity dispersion fixed~a thick 
disk velocity dispersion that increased with height would 
generate a more accurate slope. 

Modeling the disks as cylindrically rotating is only 
a crude approximation to account for the stellar cross- 
contamination. In reality, we expect the thin disk stars 
which reach large z heights to be the thin disk stars with 
larger velocity dispersions. This would mean that the 
thin disk stars at high z should also be lagging com- 
pared to the midplane thin disk stars. Ideally, we would 
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Fig. 6. — Results from fitting the midplane and offplane rotation curves as a combination of two fixed rotation curves. In the top panel, 
the rotation curves are a combination of the midplane Ha and a fiat non-rotating RC. In the bottom panel, the base rotation curves are 
the midplane Ha combined with a rotation curve counter-rotating with one-half the Ha velocity. These fits are listed in Table |4] 
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Fig. 7. — The expected thick disk lags as a function of height above the midplane and thick disk velocity dispersion. The first three panels 
show model galaxies similar to the ones in our sample. Points show the stellar lags measured from our rotation curve fits. Open points 
show lags from rotation curves where the offplane and midplane rotation curves are fit independently. Solid points show the average lag for 
the models which correct for cross-contamination of the thin and thick disk rotation curves, and are generally more reliable estimates of 
the thick disk lag. The final panel shows the results of our model when we use MW like parameters. Observed galaxies we compare to the 
models: In the upper left FGC 1642 and FGC 780; upper right FGC 1415, FGC 227, and FGC 2558; and lower left FGCE 1371. All the 
models and observations are taken at i? = 2.5/iij. FGC 1948 is excluded from the plot because there is no coherent rotation. FGC 1440 is 
excluded because we have no offplane stellar velocity measurements. 
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construct a fuUy self-consistent dynamical model of each 
galaxy, but our large uncertainties and limited LOSVD 
information would result in model degeneracies. Con- 
structing a robust self consistent dynamical model of a 
galaxy also be nefits from larger numbers of data points 
(jGirard et al.l[2006i) . With only a handful of stellar rota- 
tion curve points per galaxy, we do not have enough data 
to constrain a more complex model. We simply point out 
that when we correct for the cross-contamination of the 
rotation curves we may be over-correcting the data. We 
estimate the magnitude of the over-correction using dy- 
namical models in § [6l 

7. DUST AND PROJECTION EFFECTS 

As a final check that our observed kinematics indeed 
reflect the true stellar motions, we now explore the ex- 
pected impact of projection effects and dust extinction, 
both of which can create differences between the observed 
and underlying rotation curves. In Figure[8l we show how 
two input rotation curves are modified by being viewed 
edge-on, with and without dust. For these models, we 
assumed an exponential disk of stars and dust, and for 
simplicity only considered absorption (i.e., ignoring scat- 
tering) . The amount of dust adopted in the model would 
generate an extinction of 2.2 magnitudes in the total ap- 
parent magnitude of the galaxy. This is a rather large 
extinction for the near-IR, given that the observed galax- 
ies in our sample are only offset by 0.2 mag from the 
face-o n NIR Tully-Fisher relation (jYoachim fc DalcantonI 
|2006[). We adopted an underlying rotation curve shape 
from lCourteau (1997 ). 

As can be seen from Figure [51 the inner regions of 
the rotation curve are generally unchanged due to pro- 
jection, and the only significant changes happen in the 
outer regions, where the true rotation curve is fiat. These 
projection effects create a lag of 7.2 km s^^. 

We have not corrected our rotation curves for these 
projection effects, as we are primarily interested in the 
differences between the thin and thick disk rotation 
curves. This could lead us to make systematic errors 
in interpreting the rotation curves if the morphologies 
of the thin and thick disk are radically different, but we 
have no reason to assume this is the case. 

When dust is added to the model, it only creates an 
additional 2.6 km s^^ lag, in spite of the very high extinc- 
tion adopted here. This model is completely consistent 
with the results of lMatthews fc Wood (J2001D . who found 
that projection effects are dominant compared to extinc- 
tion in edge-on systems. We do not expect our sample 
galaxies to have larger extinctions than what is modeled 
in Figure m 

Ful l radiativ e trans f er models (iKregel &: van der KruitI 



[20051 iBianchil [2007t iXilouris et al.M 19991) . as well 



as 



comparisons of ga seous and optical rotation curves 
(|Bosma et al.lll992 ) have consistently found massive disk 
galaxies have a central face-on optical depth near unity in 
the y-band, with lower extinction levels in less massive 
syste ms like those that dominate our sample fjCalzettii 
|2001| ). Dust levels this low should not be expected to 
alter the observed rotation curve significantly, even if a 
galaxy is viewed edge-on. Moreover, most of our offplane 
rotation curves exhibit a lag compared to the midplanc. 
In contrast. If there were strong dustlanes affecting our 
midplane observations (and not the offplane), the mid- 
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Fig. 8. — Two examples of the effects dust and projection can 
have on our observed rotation curves. While projection creates 
considerable changes, the addition of dust extinction is negligible. 



plane would be the lagging component. 

The combination of working at near-IR wavelengths, 
offsetting our slit from any prominent dustlanes, and ob- 
serving intrinsically linearly rising rotation curves means 
our rotation curves should be fairly unaffected by ex- 
tinction or projection. However, the same cannot be 
said for our measured line-of-sight velocity dispersion 
(LOSVD). Unlike the rotational velocity measurement, 
which is mostly unaffected by flux contributions from 
different radii, we expect the LOSVD to be significantly 
broadened by projection effects. We also find that in 
most of our galaxies the LOSVD is very close to the in- 
strumental resolution, making any interpretation of the 
velocity dispersion suspect. Because of these challenges, 
we limit our analysis of the LOSVD to only those cases 
where we believe our measurements are of high quality 
and not dominated by the instrumental dispersion. 

8. DISCUSSION 



The results of Sections [221 and [H show that stellar kine- 
matics above the midplane display a wide range of behav- 
iors. In higher mass systems (FGCE 1371, FGCE 1498, 
FGC 1440), our midplane and offplane spectra show no 
clear signature of a hot thick disk component. The stel- 
lar rotation curves for these galaxies are well matched by 
the midplane ionized gas Ha RCs at all measured scale 
heights. All three of these galaxies converged to models 
where the rotation c urves contain no laggin g component 
(Table [H). However. I Yoachim fc DalcantonI (|2006l ) found 
that the stellar flux in higher mass galaxies are domi- 
nated by the thin disk component. Therefore, the lack 
of a significant lag in these systems is likely a result of the 
kinematically cold thin disk dominating the stellar fiux 
to scale heights of 1 kpc. This result is not completely 
unexpected, as the MW thin and thick disks should have 
simil ar luminosities 1 kpc off the midplane (jJuric et al.l 
l2008f l. We note that there is still ample photometric evi- 
dence that these higher mass galaxies contain thick disks, 
but they are simply too faint relative to the thin disk for 
modest kinematic lags to be detected spectroscopically. 

The low mass galaxies in our sample do show measur- 
able differences between the midplane and offplane ob- 
servations. At large radii, we find several galaxies where 
the offplane component is lagging compared to the mid- 
planc (Figure [7]). In three of the low mass systems (FGC 
1415, FGC 1642, and FGC 780), the lags in the offplane 
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observations become more pronounced when we correct 
for the expected thin disk contamination. These lags are 
consistent with those that are expected from dynamics 
alone (Equation [1} , provided that the thick disk has a 
radial velocity dispersion between 15 and 30 km s^^(i.c., 
10-25% of Vc). Thus, the lags in these systems do not 
necessarily require the presence of any counter-rotating 
material, although a small amount of such material could 
be present. FGC 2558 may also fall into this category; 
however the offplane RC is very similar to the midplane, 
implying this could be another galaxy where we have 
not successfully isolated the thick disk. The observed 
lags were easier to detect in these lower mass systems, 
due to their more prominent thick disks. 

The final two low mass galaxies in our sample, FGC 
227 and FGC 1948, have remarkably different rotation 
curves between the midplane and offplane. FGC 1948 
does not display coherent stellar rotation in either the 
midplane or the offplane, and therefore our subsequent 
fits converge to extreme, and probably incorrect, mod- 
els. FGC 227 does show rotation on the midplane, and a 
very low level of net rotation on the offplane. Our best 
fitting model for this galaxy has the thick disk counter- 
rotating relative to the thin disk, consistent with the 
radially increasing LOSVD which is a signature of un- 
resolved counter rotating stellar components. 

Our measurements of the LOSVD are less than enlight- 
ening. With the exception of the radial increase in the 
LOSVD in FGC 227 and the high LOSVD in FGC 1948, 
the rest of our LOSVD measurements show no signifi- 
cant trends with radius and are close to the instrumental 
resolution limit, suggesting that the radial velocity dis- 
persions of both the thin and thick disks are cold enough 
that we cannot reliably measure their velocity dispersions 
at our spectral resolution. 

Given the above results, our galaxies can be described 
as falling into three categories: The high mass systems 
which have little to no thick disk lag (or, more likely, 
thick disks which arc so faint that we have failed to 
measure their kinematics); the moderately lagging sys- 
tems; and the counter rotating system. We can now com- 
pare these results to the predictions of popular formation 
models for the thick disk. 

If thick disks are the result of gradual stochastic heat- 
ing, we would expect to always find thick disks co- 
rotating with the embedded thin disks. Moreover, with 
stronger spiral arms, larger molecular clouds, and more 
massive dark matter substructure, the high mass systems 
should be able to efficiently heat their thin disk stars into 
a thicker disk. Instead, we have found the opposite, with 
more prominent thick disks and larger lags in the lower 
mass systems, as well as evidence for counter-rotating 
stars. This seems to rule out gradual heating as the 
dominant method of thick disk formation, particularly 
for low mass galaxies. 

Forming thick disks in major mergers also does a poor 
job explaining our observations. If thick disks were pre- 
dominantly formed in major mergers that disrupt and 
heat previously thin disks, we should expect to find 
galaxies that never formed a thick disk, or that have 
failed to accrete and cool enough gas to rebuild their 
thin disk components. Major mergers also typically re- 
sult in the formation of centrally concentrated spheroidal 
components, making them a poor mechanism for forming 



thick disks in the bulgeless galaxies observed here. 

Unlike the two heating models, the variety of thick disk 
kinematics is compatible with minor mergers and/or ac- 
cretion. Presumably, the thick disk kinematics we ob- 
serve are simply the kinematics left over from the ac- 
cretion event which deposited the majority of thick disk 
stars or which triggered the formation of stars from gas 
accreted at large scale heights. The wide variety of pos- 
sible accretion events (co-rotating vs counter rotating, 
early disruption vs late disruption, high eccentricity vs 
circular initial orbit) can evolve into virialized thick disks 
with kinematics that are sometimes decoupled from the 
thin disks and that show large variation from galaxy to 
galaxy. The ubiquity of thick disks is also well explained 
by the merger/accretion scenario, given that galaxy for- 
mation in a ACDM cosmology is dominated by hierar- 
chical merging, and predicts that every galaxy has a rich 
merger history. 

Although the available data all points to a 
merger/accretion origin for the thick disk, it is difficult 
to disentangle models where thick disk stars are directly 
accreted from those where the stars form in situ further 
off th e midplane during gas rich mergers (jBrook et al.l 
[200l . 

This ambiguity results from two sources. First, there 
is no clear dividing line between what one calls a 
star-forming region off the midplane and a merging 
star- forming satellite galaxy . Second, we know from 
lYoachim fc DalcantonI (|2006[ ) that at least 75-90% of the 
baryonic accretion onto the galaxies was gaseous, and 
some fraction of this was certainly accreted in bound 
subhalos. Stars that formed initially in subhalos before 
being accreted are likely to have kinematics similar to 
those that formed from accreted gas during those same 
merging events. Presumably, one could use detailed stel- 
lar age and abundance information to help, but unfortu- 
nately this is only possible for the closest galaxies. 

There is evidence that much of the brighter inner halo 
and outer disk substructure of M31 was formed through 
accretion (iFerguson et aD2002tlKoch et al.ll2007t ). These 
features would probably resemble a thick disk if M31 
were more distant and the features were unresolved. Tak- 
ing this lesson from nearby galaxies, it is clear we are 
using smooth functions to describe thick disk that may 
actually be highly structured systems. However, the 
smooth descriptions of thick disks still provide a reason- 
able statistical description of the ensemble of accreted 
stars. 

In this study, we have measured thick disk kinemat- 
ics in only very late-type disk systems. However, thick 
disks have been photomet rically detected in a wide vari- 
ety o f Hubble tvpes (e g..[Seth et al."20G 5t iPohlen et al.l 
12004 [Morrison et al.l Il997t Ivan Dokkum et al.l I1994D . 
The kinematics in our sample are most consistent with 
merger/accretion forming for the thick disks, but, except 
for the Milky Way, there have been no measurements of 
thick disk kinematics in earlier type galaxies. 

By focusing on disk systems, we may not be sensitive 
to how thick disks form across all Hubble types. Al- 
most by definition, late-type galaxies have not suffered a 
major-merger since the formation of their stellar disks, 
otherwise they would likely possess large spheroidal com- 
ponents and be classified as an earlier type system. The 
only way pure disk galaxies could form thick disks is ei- 
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ther through accretion or stochastic heating. 

9. CONCLUSIONS 

We have expanded the kinematic observations of 
lYoachim fc DalcantonI ()2005 ') to include a total of nine 
galaxies with thick disks. Analyzing our low signal-to- 
noise spectra that contain systematic sky line residuals 
prompted us to develop a brute-force method of cross- 
correlation to extract stellar rotation curves. In galaxies 
with Vc > 120 km s"^, wc do not detect any measurable 
difference between the thin and thick disk stellar kine- 
matics. This is most likely due to a combination of thin 
disks being brighter in more massive galaxies, and the 
expected change in rotation curve as a function of scale 
height being smaller. 

In lower mass galaxies {Vc < 120 km s^^), we find a 
variety of thick disk behaviors. Thick disks are found 
with both small and large magnitude lags, including a 
counter-rotating thick disk. 

The observed kinematics are best explained by thick 
disk formation models where the thick disks in low 
mass systems are composed of stars that have been ac- 
creted from satellite galaxies or are formed at large scale 
heights from accreting gas. Models where the thick disks 
form during major mergers or through stochastic heat- 
ing seems unable to explain the wide range of thick disk 



kinematics we observe. While we strongly favor a forma- 
tion model of thick disks via accretion, we stress that this 
result can not necessarily be generalized to other Hubble 
types or higher mass systems {Vc > 120 km s~^). 
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APPENDIX 
STELLAR ROTATION CURVES IN THE PRESENCE OF SYSTEMATIC ERRORS 

Working in the ncar-IR, we find our spectra have regions which are dominated by both Gaussian and systematic 
errors caused by bright atmospheric emission lines. To properly measure stellar kinematics based on spectral absorption 
features we must employ a method that is not affected by our sky line residuals. 

There are two common techniques of deriving the kinematic information from galaxy spectra-direct y^-fitting and 
cross -correlation. In direct x^-fitting (jRix fc Whitelll992t [Kelson et al.l[2000HBarth et"aLll2002HCappellari! fc EmsellemI 
I2004D. a template star is redshifted and broad ened to fit a galaxy spectrum, while in cross-correlation techniques 
( Simkinlll974l : lTbnrv fc Davislll979tlStatleiill995D a template star is cross-correlated with the galaxy spectrum and the 
kinematic properties are deduced from the position and shape of the cross-correlation peak. 

Cross-correlation techniques have the advantage of being computationally efficient, often making use of fast Fourier 
transform algorithms. The cross-correlation technique benefits greatly from the fact that the Fourier transform of 
Gaussian noise is also Gaussian noise. In this way, noise in the galaxy spectrum transforms into random noise in the 
cross-correlation while the kinematic information becomes concentrated in a central peak. However, this is only true 
if the noise is uniform throughout the spectrum. Using a direct chi-squared fit is more computationally expensive, 
but has the added benefit of being able to weight individual wavelengths according to their specific signal-to-noise, or 
completely mask wavelengths that are affected by systematic errors. 

Although direct chi-squared fitting works well in some situations, at low S/N(<20), any direct chi-squared fitting 
routine will over-smooth the data because the low S/N continuum is best fit by a strait line (i.e., an over- broadened 
template star). In previous studies that have used direct fitting. [Kelson et al.l ([20001 ) has a median S/N of 35/A, while 
iBarth et al.l ([2002f ) report a S/N/pixel of 100-200. In contrast, our data has SNR< 20/A, due to the very low surface 
brightness of our targets. 

Because we have both low S/N and regions which require masking, we have created a fitting procedure which utilizes 
cross-correlation without making use of the computational time saving EFT techniques of previous authors. 

Traditional cross-correlation of discrete functions is defined as 



{f*9)^^J2f^ 



9t+j ■ 



(Al) 



We adopt a normalized version, where the means of the spectra have been subtracted before the cross-correlation is 
computed 



{f*9)i 






^N 



N 



(A2) 



Efc=i(/fc)^Efc=i(5fe 



where N is the number of points in the given spectra. For lags less than zero, the numerator becomes X]fe=i fk+\L\9k- 
This ensures spectra with perfectly matching shapes will have a maximum cross-correlation amplitude of unity. 
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Fig. 9. — Examples of cross-correlating in the presence of different types of noise. In the left hand column, we show a model galaxy 
spectrum (top) and stellar template (bottom). In the right hand column, we plot the galaxy-star cross-correlation (solid) and stellar 
auto-correlation (dotted) and note the velocity error resulting from comparing the two. (a) The ideal case of a high signal-to-noise galaxy 
spectrum, (b) Results from a galaxy spectrum with a S/N/A~ 10. (c) Spectra with a small region of very low S/N affecting a section of 
one of the Ca absorption features, similar to how bright sky lines leave residuals on our spectra, (d) A traditional cross-correlation where 
the noisy region has been set to the continuum, (e) Our new cross-correlation technique where we compute the cross-correlation excluding 

* "PHISn^^ 'i!v^°flefine masks S for each spectrum which have values of 1 in regions of good data and for masked 
wavelengths. Given a stellar spectrum S and Galaxy spectrum G that are binned in logarithmic wavelength intervals 
and have both been normalized by division of a low order polynomial and had their means subtracted, we compute 
our modified cross-correlation as 



iS*G)i 
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We then generate a model galaxy spectrum M by redshifting and broadening the stellar template, M{x) = S{x + 
v) (X) B{x) where B{x) is a Gaussian broadening function, w is a velocity shift, and (X> represents convolution. We then 
calculate the model's modified cross-correlation using the masks from the actual galaxy spectrum 



^N 



N 



EtiiSkSf^soy^k 



,{M,5p<^) 



We vary the velocity shift and broadening to minimize the x^ between {S -kG) and (S-kM). We focus on the region 
of the primary peak, and clip regions beyond the bracketing local minima. Examples of traditional cross-correlation 
and our modified cross-correlation are shown in Figure [9l In general, our masked cross-correlation technique cannot 
reproduce the excellent fits that are possible with data that is unaffected by systematics, but we can reduce the errors 
to be of order 5 km s^^ in our typical spectra. 
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